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ABSTRACT 

In this paper we create a complete gravity torque map of the disk of the LINER/Seyfert 1.9 galaxy NGC 4579. We quantify the 
efficiency of angular momentum transport and search for signatures of secular evolution in the fueling process from r ~ 15 kpc down 
to the inner r~50pc around the Active Galactic Nucleus (AGN). We use both the 1-0 and 2-1 line maps of 12 CO obtained with the 
Plateau de Bure Interferometer (PdBI) as part of the NUclei o/Gato/^-(NUGA)-project. The CO(l-O) and CO(2-l) PdBI maps, at 
2.0" x 1.3" and 1.0"x0.6" resolution respectively, include short spacing correction provided by IRAM-30m data. We derive the stellar 
potential from a NIR (K band) wide field image of the galaxy. The ^T-band image, which reveals a stellar bar, together with a high 
resolution HI map of NGC 4579 obtained with the Very Large Array (VLA), allow us to extend the gravity torque analysis to the 
outer disk. Most of the molecular gas mass in the inner r <2 kpc disk is distributed in two spiral arms, an outer ring, and a central 
lopsided disk. The bulk of the gas response traced by the CO PdBI maps follows the expected gas flow pattern induced by the bar 
potential in the presence of two Inner Lindblad Resonances (ILR) at r ~500 pc and r ~1.3 kpc . We also detect an oval distortion in 
the inner r ~200 pc of the AT-band image. The oval is not aligned with the large-scale bar, a signature of dynamical decoupling. The 
morphology of the outer disk suggests that the neutral gas is currently piling up in a pseudo-ring formed by two winding spiral arms 
that are morphologically decoupled from the bar structure. The pseudo-ring is located inside the bar corotation (t C r ~6 kpc) and close 
to the predicted position of the Ultra Harmonic Resonance (UHR) (r UHR ~3.8 kpc). 

The derived gravity torque budget in NGC 4579 shows that the fueling process is at work on different spatial scales in the disk. In 
the outer disk (r >2 kpc), the decoupling of the spiral allows the gas to efficiently populate the UHR region, and thus produce net 
gas inflow on intermediate scales. The corotation barrier seems to be overcome due to secular evolution processes. The gas in the 
inner disk (r <2 kpc) is efficiently funneled by gravity torques down to r ~300 pc. Closer to the AGN (r <200 pc), gas feels negative 
torques due to the combined action of the large-scale bar and the inner oval. The two m - 2 modes act in concert to produce net 
gas inflow down to r ~50 pc, providing a clear smoking gun evidence of fueling with associated short dynamical time-scales (~l-3 
rotation periods). 

Key words. Galaxies individual: NGC 4579 - Galaxies:ISM - Galaxies: kinematics and dynamics - Galaxies: nuclei - 
Galaxies :Seyfert - Radio lines: galaxies 



1. Introduction 

Nuclear activity in galaxies is understood as the result of the 
feeding of supermassive black holes (SMBH). SMBH are sus- 
pected to be a common component in most, perhaps all, galaxies 
with a significantly massive bulge (e.g., see review by Ferrarese 
& Ford 2005 and references therein). Active Galactic Nuclei 
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(AGN) must be fed with material coming from the disk of 
the host galaxy. The supplying gas, lying originally far away 
from the gravitational influence of the black hole, must lose 
virtually all of its angular momentum in the fueling process. 
While for high luminosity AGNs, a good correlation between 
the presence of ~kpc-scale non-axisymmetric perturbations and 
the existence of activity has been found (e.g., Hutchings & 
Neff I1992L the case for a similar correlation in low luminosity 
AGNs (LLAG Ns) is weak if any (e.g ., Mul chaey & ReganH997l 
Knapen et al. I20UU1 Marquez et al. l2000l Combes [2003]). The 
search for a universal mechanism for AGN feeding in LLAGNs 
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is probably complicated by the fact that the AGN duty cycle 
(~10 6-7 years) might be shorter than the lifetime of the feeding 
mechanism itself (e.g., Wada 2004; Hopkins & Hernquist 2006; 
King & Pringle 2007 ). In spite of all the theoretical and obser- 
vational efforts, finding a solution for the LLAGN fueling prob- 
lem has thus far remained an elusive aim (e.g., see review by 
Martini |3D@B. 

The study of the content, distribution and kinematics of in- 
terstellar gas is a key to understand the origin and maintenance 
of nuclear activity in galaxies. The processes involved in AGN 
fueling encompass a wide range of scales, both spatial and tem- 
poral, which have to be studied. Probing the gas flow from the 
outer disk down to the central engine of an AGN host, requires 
the use of specific tracers of the interstellar medium adapted to 
follow the change of phase of the gas as a function of radius. 
As most of the neutral gas in galactic nuclei is in the molecu- 
lar phase, low-J rotational lines of carbon monoxide (CO) are 
the best choice to undertake high-spatial resolution (-1") inter- 
ferometer mapping of the central kiloparsec disks of AGNs. On 
the other hand, the 21cm line emission of atomic hydrogen (HI) 
is best suited to trace the distribution and kinematics of neu- 
tral gas in the outer disk reservoirs of AGNs (r >5 kpc) with 
moderate spatial resolution (^10-20") using an interferometer. 
Besides providing a handle on the content and distribution of the 
gas, CO and HI interferometer maps give a sharp view of the gas 
kinematics. The combined information extracted from CO and 
HI is essential to characterize gravitational instabilities and con- 
strain theoretical models of the gas flows at all spatial scales in 
galaxy disks. 

The NUclei of GAlaxies-NUGA-project, described by 
Garcfa-Burillo et al. (2003a, 2003b), is the first high-resolution 
(~0.5"-l") 12 CO survey of 12 nearby (D=4-40 Mpc) LLAGNs 
including the full sequence of activity types (Seyferts, LINERs 
and transition objects). Observations, carried out with the IRAM 
Plateau de Bure Interferometer (PdBI), have been completed 
early 2004. NUGA surpasses in both spatial resolution (10- 
100 pc) and sensitivity (3<x-detection limit^a few 10 5_6 Mq) the 
previous surveys of nearby AGN conducted at the Owens Valley 
Radio Observatory-OVRO (Baker l2000l Jogee et al. l200TT) and 
at the Nobeyama Radio Observatory-NRO (Kohno et al. 12001b . 
NUGA maps allow us to probe the gas flows at critical spa- 
tial scales (<100 pc) where secondary modes embedded in the 
kpc- scale perturbations are expected to take over in the fueling 
process. A bottom line result of NUGA is the identification of 
a wide range of morphologies in the central kpc-disks of these 
LLAGNs. This includes one-arm spirals or m - 1 instabilities 
(NGC 4826: Garcfa-Burillo et al. I2003bi NGC 3718: K rips et 
al. 120051) . symmetric rings (NGC 7217: Combes et al. 120041 
NGC 3147: Casasola et al. 2007), as well as gas bars and two- 
arm spirals (NGC 4569: Boone et al. 120071 NG C 2782: Hunt 
et al. 12001 NGC 6574: Lindt-Krieg et al. 120081) . As such, this 
result is suggestive of an evolutionary scenario where several 
mechanisms cooperate to feed the central engines of LLAGNs. 
We have performed a detailed case-by-case study of the distri- 
bution and kinematics of molecular gas in the galaxies of our 
sample, and interpreted these in terms of evidence of ongoing 
feeding. The gas response to the stellar potential is characterized 
with the help of high-resolution optical and NIR images of the 
galaxies. We adopt this in-depth approach to take full advantage 
of the high quality of NUGA maps. 

In a previous pilot study, Garcfa-Burillo et al. (12005b (here- 
after GB05) analyzed the efficiency of stellar gravity torques 
to drain the gas angular momentum in the central kiloparsec 
of four NUGA targets: NGC 4321, NGC 4826, NGC 4579 and 



Table 1. Observational parameters of NGC 4579. 



Parameter 


Value 


Reference 


<^/2ooo (dynamical center) 


12 /z 37 m 43.52 5 


This work 


S 72000 (dynamical center) 


11°49'05.5" 


This work 




1466 km s _1 


This work 


RC3 Type 


SAB(rs)b 


NED 


AGN Type 


LINER/Syl.9 


NED 


t Type 


2.8 


LEDA 


Inclination 


36° 


GB05 


Position Angle 


95° 


GB05 


Distance 


19.8 Mpc (1"=97 pc) 


GB05 


M B 


-21.68 mag 


LEDA 


M(HI) 


1.7 xlO 8 M Q 


Haan et al. (I2008al> 


M(H 2 ) 


5 xlO 8 M Q 


This work 


Lfir 


9.6 xlO 9 L Q 


Sanders et al. (2003 ) 



NGC 695 1 . Results of this analysis indicate that, paradoxically, 
feeding due to the stellar potential is presently quenched close 
to the AGNs. Due to the inhibiting action of bars (NGC 4321, 
NGC 4579) or weak oval perturbations (NGC 4826, NGC 695 1), 
the derived gravity torque budget is seen to be positive inside 
r ~ 150-200 pc in the four LLAGNs analyzed in this paper. This 
situation is consistent with the theoretical picture where gas gets 
trapped in the Inner Lindblad Resonance (ILR) of large-scale 
bars. ILRs favor the feeding of a starburst in nuclear rings but can 
halt the gas flow inwards (Combes fT988l Regan & Teuben [2004b . 
GB05 speculate that the agent responsible for driving inflow 
at smaller radii may be transient (lifetimes 10 6 yrs); finding a 
smoking gun evidence of fueling would be thus difficult (see 
however the case of NGC 2782 analyzed by Hunt et al. 12505b . 
As an alternative explanation, GB05 estimate on a case-by-case 
basis that the gravity torque barrier associated with the ILRs in 
these galaxies could be overcome by viscosity. In this scenario, 
gravity torques and viscosity coordinate efforts to produce recur- 
rent episodes of activity during the typical lifetime of any galaxy. 
The weakening of bars, triggered by the radial re-distribution of 
gas and the ensuing angular momentum exchanges (Bournaud & 
Combes [2002b . may allow AGN feeding episodes to occur. This 
would account for the presence of filamentary dusty structure 
connecting nuclear rings to the central engines of some LLAGNs 
(e.g., Peeples & Martini 12006b . 

NGC 4579 is a SAB(rs)b galaxy classified as an intermediate 
type 1 object (LINER/Seyfert 1.9) by Ho et al. (TT997b . The AGN 
nature of NGC 4579 is confirmed by the detection of broad H a 
and UV lines (Sta ufferggH Filippenko & Sargent[[985j Barth 
et al. 119961 I2U0TI Maoz et al. [TW8l It also has an unresolved 
nuclear hard X-ray (variable) source with a prominent broad Fe 
K a line (Terashima et al. 120001 Ho et al. 1200 lai Eracleous et 
al. 120021 Dewangan et al. 12004b . A non-thermal radio contin- 
uum source is detected at the position of the AGN (Hummel et 
al.H987; Ho & Ulvestad|2003B Ulvestad & Ho l20OT1: Krips et 
al. [2007b . Hubble Space Telescope (HST) optical and UV im- 
ages resolve the central kiloparsec disk of NGC 4579 and show 
a nuclear component surrounded by a highly structured disk 
(Pogge et al. 120001 Contini |2004b . The galaxy is part of the Virgo 
cluster, although neither the morphology nor the kinematics of 
the HI gas disk, mapped by Haan et al. (2008a), show any signs 
of significant stripping at work. 

In this paper we re-determine the gravitational torque budget 
in the central kiloparsec of NGC 4579 using both the 1-0 and 
2-1 line maps of 12 CO obtained inside the NUGA project. The 
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Fig. 1. a) (L^/r panel) CO(l-O) line emission profiles observed with the IRAM 30m telescope in NGC 4579. We have observed with 
a 5x5-point grid and a 7" spacing. Offsets are relative to the AGN. Velocity scale goes from -400 km s" 1 to 400 km s" 1 with respect 
to u=t;o=1520 km s" 1 in LSR scale, and temperature scale (in T*) goes from -30 to 100 mK. b) (Right panel) Same as a) but for the 
CO(2-l) line. Same velocity scale as in a); temperature scale from -40 to 150 mK. All spectra have been smoothed to a resolution 
of 10 km s" 1 . 



new maps used in this work include the short spacing correc- 
tion provided by the IRAM- 30m data obtained in the two lines 
of CO, and thus give a more reliable estimate of the implied gas 
flow time-scales, compared to the first estimate of GB05. We de- 
rive the stellar potential from a NIR (K band) wide field image 
of the galaxy disk; this allows us to minimize the residual ef- 
fects of extinction that were present in the HST /-band image 
that we originally used in the first gravity torque analysis of this 
galaxy (GB05). Furthermore, the large field-of-view of the K- 
band image, together with the availability of a high resolution 
and sensitivity HI map of NGC 4579, recently obtained by Haan 
et al. (I2008ab . allow us to extend the gravity torque analysis to 
the outer disk of the galaxy. A complete gravity torque map of 
NGC 4579, makes it possible to quantify the efficiency of angu- 
lar momentum transport from r ~ a few kpc down to the inner 
r~50pc around the AGN and search for signatures of secular 
evolution in the accretion process. This is a key to test the fuel- 
ing scenario described by GB05. 

We will assume a distance to NGC 4579 of around 20 Mpc 
for consistency with Garcfa-Burillo et al. (12005b . which im- 
plies 1" ^lOOpc. This is a good compromise between the differ- 
ent values reported in the literature ranging from Z)=16.8 Mpc 
(Tully & Fisher [1988]) to £>=22.6 Mpc (LEDA, 2006). Similarly 
to GB05 we will assume that the inclination angle and the po- 
sition angle of NGC4579's disk are /=36° and PA=95°, respec- 
tively, in rough agreement with previous determinations (e.g., 
Rubin et al. [T9991 Koopmann et al. I200T1 Daigle et al. [2D06b . 
The basic observational parameters of NGC 4579 are listed in 
Table[Q 

We describe in Sect.[2]the observations used, including high- 
resolution CO and HI maps as well as ground-based NIR images 
of NGC 4579. Sects. [3] to [6] describe the distribution and kine- 
matics of molecular gas. We describe the stellar structure and 
star formation of NGC 4579 in Sect. [TJ Sect. [8] computes from 
NIR images the gravitational potentials and forces, and deduces 



from the combined CO and HI maps the effective torques applied 
to the gas. From these torques, it is possible to derive time-scales 
for gas flows and discuss whether gravity torques alone are effi- 
cient enough to feed the current activity in NGC 4579. The gen- 
eral implications of these results for the current understanding of 
AGN feeding are summarized in Sect. [£l 

2. Observations 

2.1. Interferometer CO maps 

Observations of NGC 4579 were carried out as part of the 
NUGA survey conducted with the PdBI during two years and 
finished by February 2003. We used the ABCD set of config- 
urations and the six antennas of the array in dual frequency 
mode to reach the highest spatial resolution (<1" at the high- 
est frequency) but also to maximize sensitivity to all spatial fre- 
quencies in the maps. We observed simultaneously the emission 
of the J=l-0 and J=2-l lines of CO in single fields of sizes 
42" and 21", respectively, centered at a /2 ooo=12 / *37 m 43.58* and 
£/2000=ll°49'02.5". The first results derived from the analysis 
of the CO(2-l) maps were discussed by GB05. During the ob- 
servations the spectral correlator was split in two halves cen- 
tered at the transition rest frequencies corrected for the assumed 
recession velocity v (LSR)= 1520 km s" 1 . The correlator con- 
figuration covers a bandwidth of 580 MHz for each line, using 
four 160 MHz- wide units with an overlap of 20 MHz; this is 
equivalent to 1510km s" 1 (755km s" 1 ) at 115GHz (230GHz). 
Visibilities were obtained using on- source integration times of 
20 minutes framed by short (~ 2 min) phase and amplitude cal- 
ibrations on the nearby quasars 1156+295 and 3C273. The ab- 
solute flux scale in our maps was derived to a 10% accuracy 
based on the observations of primary calibrators whose fluxes 
were determined from a combined set of measurements obtained 
at the 30m telescope and the PdBI array. The bandpass calibra- 
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Fig. 2. 12 CO velocity-channel maps observed with the PdBI in the nucleus of NGC 4579 with a spatial resolution of 2.0 // xl.3 // at 
PA=26 ° (beam is plotted as a filled ellipse in the bottom left corner of each panel). We show a field of view of 42", i.e. -the diameter 
of the primary beam at 115 GHz. The phase tracking center is indicated by a cross at #2000= 12 / *37 m 43. 58* and £2000=1 1°49 / 02.5 // . 
Velocity-channels are displayed from v =200 km s" 1 to v =-260 km s" 1 in steps of 10 km s" 1 . Velocities are in LSR scale and relative 
to f =l»o=1520 km s" 1 . Contour levels are -3<x, 3<x to 24<x in steps of 3<x where the 1-sigma rms cr=2.3 mJy beam -1 . There are very 
few pixels with negative fluxes below -3cr inside the field-of-view. 



tion was carried out using 3C273 and is accurate to better than 
5%. The point source sensitivities derived from emission-free 
channels of 10 km s" 1 width are 2.55 mJy beam -1 in CO(l- 
0) and 5.44 mJy beam -1 in CO(2-l). Images of the continuum 
emission of the galaxy at 115 GHz and 230 GHz have been ob- 
tained by averaging those channels free of line emission. The 
image reconstruction was done using standard IRAM/GAG soft- 
ware (Guilloteau & Lucas 2000 ). We used natural weighting and 
no taper to generate the 1-0 line maps with a field of view of 
76.8" and 0.15" sampling; the corresponding synthesized beam 
is 2.0" x 1.3", PA=26°. We also used natural weighting to gener- 
ate 2-1 maps with a field of view of 51.2" and 0.10" sampling; 
this enables us to achieve a spatial resolution <\" (1.0" x 0.6", 
PA =25°). The conversion factors between Jybeam -1 and K are 
36 K Jy" 1 beam at 1 15 GHz, and 40 K Jy" 1 beam at 230 GHz. 



2.2. IRAM-30m CO observations and short spacing 
correction 

With the aim of including the short spacing correction in the fi- 
nal image, we have mapped the emission of the 2-1 and 1-0 
lines of CO in NGC 4579 with the IRAM 30m telescope. We 
have covered the central 28" (2.8 kpc) region of the galaxy us- 
ing a grid of 5x5 points with a 7" spacing. The area covered 
is larger than the primary beam size of the PdBI at 230 GHz, 
and a significant fraction (-70% ) of the corresponding field of 
view at 115 GHz. Observations were carried out in two observ- 
ing runs between July 2002 and June 2004. At 115 GHz and 
230 GHz, the telescope half-power beam widths are 22" and 
12", respectively. The 3 mm and 1 mm receivers of the 30m tele- 
scope were tuned to the redshifted frequencies of the lines. The 
backends were two 1 MHz filter banks and auto-correlator spec- 



S. Garcfa-Burillo et al.: Molecular gas in NUclei of GAlaxies (NUGA) 




■f 




.+ 



+ 



+> 



■ ,+ 



176 ' 



+ 



' I H P I | I 



0_° 



.+ 



•+ 



■ I l i I 



A, 



12 



CO(2-l) 



12 n 37 m 44 s Oi3 s 5 43 s 
R.A. 



Fig. 3. Same as Fig. 1 but for the 2-1 line of CO. Spatial resolution reaches 1.0"x0.6" at PA =-155° (beam is plotted as a filled 
ellipse in the bottom left corner of each panel). We show a field of view of 21", i.e. -diameter of the primary beam at 230 GHz. 
Velocity-channels are displayed from v=200 km s" 1 to v=-260 km s" 1 in steps of 10 km s" 1 , with same reference as used in Fig. 
1. Contour levels are -3<x, 3<x to 19<x in steps of 2<x where the 1-sigma rms <x=4.6 mJy beam" 1 . There are very few pixels with 
negative fluxes below -3<x inside the field-of-view. 



trometers. The velocity range covered was 1340 km s" 1 for the 
3 mm lines and 670 - 1200 km s" 1 for the 1 mm lines. Typical 
system temperatures during the observations were -250 K at 
3 mm and -600 K at 1 mm. All receivers were used in single 
side-band mode (SSB), with a high rejection of the image band: 
>15dB at 1mm and >20dB at 3 mm. The latter assures that 
the calibration accuracy for the bulk of our data is better than 
20%. Throughout the paper, velocity-integrated line intensities 
(7) are given in antenna temperature scale, T*. The r* scale 
relates to the main beam temperature scale, T m b, by the equa- 
tion r m b = (F Q ff/B Q ff)T*, where F q q and # e ff are the forward 
and beam efficiencies of the telescope at a given frequency. For 
the IRAM 30m telescope F eff /B eff = 1.27 (1.75) at 115 GHz 
(230 GHz) and S/T mh = 4.95 JyK" 1 . Wobbler switching mode 
was used to assure flat baselines, taking a reference position off- 
set by 4 r in azimuth. The pointing was regularly checked on con- 
tinuum sources and the accuracy was 3" rms. 



Short spacings were included in the combined PdBI+30m 
image using the SHORT-SPA CE task available in the GILDAS 
software (Guilloteau & Lucas 120001) . The relative weights are 
chosen so as to guarantee that similar absolute weights are taken 
for the single-dish data and the interferometer data within a ring 
in the UV plane going from 1.25 D/A to 2.5 D/A(D= 15 m). In 
the combined map we ended up recovering all the missing flux 
and simultaneously kept the spatial resolutions very close to 
those of the original PdBI maps. Throughout this paper we use 
the combined PdBI+30m data cube, i.e., corrected by short spac- 
ings. We estimate that the map including only the PdBI spac- 
ings, published by GB05 recovers -65% of the total CO(1-0) 
flux measured with the combined PdBI+30m map. A similar re- 
trieval factor is derived for the CO(2-l) line. The 1-cr noise lev- 
els in 10 km s _1 -wide channels for the PdBI+30m images are 
2.3 mJy beam" 1 in CO(1-0) and 4.6 mJy beam" 1 in CO(2-l). 
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Fig. 4. a)(Left panel) The continuum map at 1 mm obtained with the PdBI (in contours from 20% to 100%, in step of 10% of the max- 
imum=10 mJy beam -1 ) is overlaid on the VLA radio-continuum image at 6 cm (in grey scale) taken from Ho & Ulvestad ( 2001b). 
The AGN is clearly identified with a non-thermal point source. (Aa, A£)-offsets are here relative to the phase tracking center. 
b)(Middle panel) We show the emission associated with the AGN source (torus) observed in the CO(1-0) line. c)(Right panel) The 
same as b) but for CO(2-l), centered around vq. The derived systemic velocity v sys =v-vo=-50 km s _1 . 



The 30m maps probe the CO emission of molecular gas in 
the central r~1.4 kpc of NGC4579 (Fig. [T]). The emission is 
spatially resolved in CO(2-l) and is suggestive of an elongated 
gas distribution, mostly along the east- west axis. This picture 
is confirmed by the higher resolution interferometer data (see 
Sect. 11). 



2.3. HI maps 



The 21cm HI line emission in NGC 4579 was observed in 2003 
and 2004 using the NRAO Very Large Array (VLA) in its C 
and D-configurations. A first analysis of the results obtained in 
NGC 4579 have been published by Haan et al. (I2008al) as part 
of a comprehensive spectroscopic imaging HI survey conducted 
in a sample of 16 LLAGNs including all NUGA galaxies. This 
HI survey reaches a moderate spatial resolution (-30-40") and 
high sensitivity (3-cr detection limit of <10 19 cm -2 ). Calibration 
procedures followed are fully described in Haan et al. (2008a). 
The data issued from this HI survey have been recently the basis 
of an overall study of the in/out-flow rates in 7 NUGA galax- 
ies (Haan et al. I2008bb . The HI maps of NGC 4579 used in 
this paper have been produced by robustly weighted imaging 
of the data in order to maximize spatial resolution while keep- 
ing a high sensitivity threshold. With this procedure we reach a 
spatial resolution of 23.1" x 21.6" and a 3-cr sensitivity limit 
of 0.51 mJy beam" 1 . Analysis was done with the Groningen 
Image Processing SYstem (GIPSY) as described by Haan et 
al. d200 8a). The channel maps were combined to derive an HI 
intensity map with a flux cut-off of 3<x. The HI emission in the 
disk of NGC 4579 lies well inside the primary beam field of view 
of 30' , therefore no attempt has been made to correct for primary 
beam attenuation. 



2.4. Optical, UV and IR images 

We used several optical, UV and IR images of NGC 4579 to 
study the star formation pattern and the stellar structure of the 
disk of the galaxy as listed below. 

We first acquired from the HST archived broadband im- 
ages of NGC 4579, including two WFPC2 images (F547M and 
F791W). Similarly, a continuum- subtracted H a emission line 
image was created by subtracting the associated continuum-band 
image, as described in Pogge et al. (2000). The optical images 
were combined using (crreject) to eliminate cosmic rays, and 
calibrated according to Holtzman et al. (1995). The "pedestal" 
effect (see Boker et al. 119991) was removed with the van der 
Marel algorithm^. Sky values were assumed to be zero since the 
galaxy filled the WFPC2 frames, an assumption which makes an 
error of ~ 0.1 mag at most, in the corner of the images. 

We used the ultraviolet (UV) image of NGC 4579 obtained 
with the Hubble Space Telescope's Advance Camera for Surveys 
(ACS) and its detector High Resolution Channel (HRC) that 
provide -0.028" x 0.025" /pixel spatial resolution. This image 
correspond to the F330W band whose central wavelength is 
~ 3300 A. (Maoz et al. 2 005 ). The image is available already 
reduced with the Space Telescope Science Institute (STScI) 
pipeline. 

We also used the far-ultraviolet (FUV) image of the GALEX 
satellite, whose band is centered at A e ff = 1516 A. This image 
has been obtained with a total exposure time of 1586 sec and 
covers a square region on the sky of size ~ 5800" x 5800". The 
image is available already reduced with a slightly modified ver- 
sion of the GALEX data pipeline (Gil de Paz et al. 12007b , ex- 
pressed in intensity units and sky-subtracted. 

In order to derive the stellar potentials, we adopted the K- 
band image obtained using the INGRID NIR camera of the 

1 Based on observations made with the NASA/ESA Hubble Space 
Telescope, obtained from the data archive at the Space Telescope 
Science Institute. STScI is operated by the Association of Universities 
for Research in Astronomy, Inc. under NASA contract NAS 5-26555. 

2 http://www. stsci. eduf marel/ 'software/pedestal html 
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Fig. 5. a) The CO(l-O) integrated intensity map obtained with 
the PdBI in the nucleus of NGC4579. The map is shown in 
grey scale with contour levels 0.2, 0.5, 1, 1.7, 2.6, 3.7 and 
5 Jy km s" 1 beam -1 . The filled ellipse at the bottom left cor- 
ner represents the CO beam size. (Aar, A£)-offsets are with re- 
spect to the location of the AGN (marked by the star): (RA2000, 
Dec 2 ooo)=(12 / *37 m 43.52*, 11°49'05.5"). b) Same as a) but here 
for the CO(2-l). Contours levels from 0.3, 0.6, 1 to 9 in steps 
of 1 Jy km s" 1 beam -1 . Positive intensity contours have been de- 
rived using a 3-cr clipping on both data cubes. The position of the 
N-arm, S-arm and O-ring, as defined in the text, are highlighted. 



William Herschel Telescope (WHT) by Knapen et al. (2003). 
The field of view provided by the camera is 4.2' x 4.2' with a 
pixel size of 0.242". The details of the data reduction and cal- 
ibration procedures for this image are described in Knapen et 
al. (120031) and references therein. 

IRAC images in bands 3.6 jim and 8 jjm were retrieved from 
the Spitzer archive of NGC 4579 (Kennicutt et al. 120031 Theses 
images cover a region of the sky of size ~ 1400" x 1500" and have 
pixel size of 0.75". Images are reduced with the SINGS IRAC 
pipeline and are calibrated in MJy sr" 1 units. We have recentered 
on the AGN and deprojected onto the galaxy plane both images 
to superpose these IR images to the HI map. 




Aa(arcsec) 

Fig. 6. a) The CO(1-0) integrated intensity map obtained with 
the PdBI (contour levels are 0.2, 0.5, 1.0, 1.7, 2.6, 3.7 and 
5Jy km s" 1 beam -1 ) is overlaid on the V - I color image 
from HST (color scale) observed in the nucleus of NGC 4579. 
(Aa, A£)-offsets are with respect to the location of the AGN 
(marked by the star), b) Same as a) for the CO(2-l) line 
with contour levels from 0.3, 0.6, 1.0 to 9 in steps of 
1 .0 Jy km s" 1 beam -1 . A close-up view of the inner 200 pc region 
shows a lopsided disk with the AGN lying on its southwestern 
edge. 



3. Dynamical center and position of the AGN 

Figures |2]and[3]show the velocity-channel maps of CO(1-0) and 
CO(2-l) emission in the central region of NGC 4579, respec- 
tively. The observed kinematical pattern is characteristic of a 
spatially resolved rotating circumnuclear disk of ~2 kpc-radius. 
The gas kinematics indicate that the eastern (western) side of 
the CO disk is red (blue)-shifted with respect to the reference 
velocity v initially assumed in this work. This is in agreement 
with the kinematical major axis of the circumnuclear disk be- 
ing oriented roughly along the RA-axis, an orientation similar 
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to that of the outer disk (r>5-10 kpc) probed by HI and H a 
emission (Guhathakurta et al. |1988t Chemin et al. 2006; Haan et 
al. 12008 al) . The global kinematics are perturbed by non-circular 
motions (see Sect. [6]). However, and contrary to previous claims 
(Chemin et al. 2006 ), we do not find evidence that the kinemati- 
cal major axis of the inner disk and that of the outer disk are sub- 
stantially different. In this paper we adopt a common PA =95° at 
all radii in rough agreement with the values derived from H a and 
HI (see Sect. ED. 

The best fits for the center of symmetry of rotation 
and systemic velocity are (Aar,AJ)=(-0.9", 3") and v^f=v Q - 
50 km s _1 = 1470 km s" 1 , respectively. Within the errors, the de- 
rived CO dynamical center (a/2000 = 12 h 37 m 43.52 s , £72000 = 
ll°49 / 05.54 // ) coincides with the position of the AGN, as de- 
termined by multi-frequency very long based interferometry 
(VLBI) observations of the non-thermal continuum emission of 
NGC4579 (Anderson et al. l2004l Krips et al. 120571 . This coin- 
cidence is confirmed by our detection of a point source in con- 
tinuum emission at 115 GHz (with a flux of 11.5+0.3 mJy) and 
230 GHz (with a flux of 11.0+1 mJy) (see Fig.0). As fully dis- 
cussed by Krips et al. (2007 ), the spectral index of the continuum 
emission from 1 GHz to 230 GHz indicates that the spectrum of 
the AGN source is flat or slightly inverted, as frequently found 
inLLAGNs. 

Weak but still statistically significant emission has also been 
detected in the 2-1 and 1-0 lines of CO coming from a spa- 
tially unresolved source at the position of the AGN. Emission 
is detected at a >4-<x level in integrated intensities over a ve- 
locity range Ay -150 km s" 1 . This suggests the existence of 
a molecular torus/disk of radius r<40pc, an upper limit deter- 
mined by the beam size at 230 GHz. The systemic velocity de- 
rived from CO(y^=1470±10 km s" 1 =^f L =1466+10 km s" 1 ) 
is -50 km s" 1 blueshifted with respect to the value determined 
from HI (^f L =1519+10 km s" 1 from RC3). The value of v sys 
based on CO, besides symmetrizing the velocity field of the 
gas inside r ~2 kpc, coincides with the velocity centroid of the 
CO(2-l) line emission detected in the torus/disk. The 1-0 spec- 
trum towards the AGN is slightly blueshifted with respect to v sys 
due to beam smearing effects at this frequency. 

4. Molecular gas distribution 

The velocity-integrated CO(1-0) flux obtained in the combined 
PdBI+30m map within the 42" primary beam field of the PdBI 
is S co = 150 Jy km s" 1 . T his is -90% of the total flux detected by 
Kenney & Young ( 1989 ) towards the central position of a raster 
map of NGC 4579, made along the major axis of the galaxy 
with the Five College Radio Astronomy Observatory (FCRAO) 
telescopeB Assuming a CO-to-H2 conversion factor typical 
of galaxy nuclei X= N(H 2 ) /I C o=23xl0 20 cm" 2 K" 1 km" 1 s 
(Solomon & Barrett 11991b , the total H2 mass derived from 
the PdBI+30m map within the 42" PdBI field of view 
is M(H 2 ) ~5xl0 8 M©. Including the mass of helium, the 
corresponding total molecular gas mass is M gas =M(H 2 + 
He)=136xM(H 2 ) ~7xl0 8 M o . The distribution of molecular 
gas is best seen in Figs.[5]2,Z?, which show the CO intensity maps 
obtained by integrating the emission in velocity channels from 
v-v sys = -210 to 260 km s" 1 using a 3-cr clipping. 

The sensitivity of the CO(1-0) map is reduced beyond the 
primary beam radius at this frequency r =21" (~2 kpc). This ex- 

3 Note that the two flux estimates are directly comparable as the 
FWHM of the FCRAO at 115 GHz is 45", i.e., similar to the corre- 
sponding PdBI primary beam size at this frequency 



plains why we do not see in our maps the large-scale two-arm 
spiral structure detected in the single-dish CO(1-0) map of Kuno 
et al. (12007b . This large-scale spiral develops from r ~4 kpc to 
r ~ 10 kpc and shows a good correspondence with the HI pseudo- 
ring detected by Haan et al. (2008a) (see Fig. [Tib). Normalized 
to a distance D=20 Mpc, we estimate that the total molecular gas 
mass inside r -10 kpc amounts to M gas ~4x10 9 Mq. 

Inside r~2 kpc we identify in the PdBI maps three main 
components in the distribution of the 2-1 and 1-0 line emission: 

- The central disk: the 2-1 line map of Fig. 0? gives a sharp 
view of the distribution of molecular gas at radii r < 200 pc. 
As shown in Figs. [6] and [TJ there is little molecular gas at 
r<200 pc from the central engine. The most massive gas 
complexes, of a few- 10 6 M o , lie east and north of the nucleus 
at r~150pc (complexes E and N in Fig|7j respectively). 
Complexes E and N are connected by an arc of low-level 
CO(2-l) emission (complexes E' and N' in Fig. [7]). These 
constitute altogether the molecular gas counterpart of a lop- 
sided disk that can be clearly identified as a dusty feature in 
the V-I color HST image of the galaxy (Fig. [6]). The AGN 
lies close to the southwestern edge of the central disk. As 
discussed in Sect. [3j there is also weak CO emission related 
to the AGN; this corresponds to a molecular gas torus or disk 
of M gas =M(H 2 + He) ~10 6 M o . 

- The inner spiral arms: the bulk (>80%) of the molecular 
gas as traced by the two CO lines is piled up in two highly 
contrasted spiral lanes that lie at the leading edges of the 
~12 kpc-diameter bar, identified in the ^-band image of 
Knapen et al. (2003)(see also Eskridge et al. 2002 and Jarrett 
et al. 2003). The described geometry is also illustrated in 
Fig. [Til?, c. which shows the deprojected version of the CO 
maps overlaid on the ^T-band image of the galaxy. The CO 
arms develop from r ~300 pc to r ~1.3 kpc. At closer exami- 
nation, it can be seen that the CO spiral lane stretching north 
of the AGN (hereafter denoted as N-arm) is more regular 
and better delineated than its southern counterpart (hereafter 
referred to as S-arm). Moreover, the inner segment of the N- 
arm is significantly closer to the AGN compared to that of the 
S-arm (Fig. [6]). As discussed in Sect. [7] (see also Fig. [12]), a 
similar asymmetry is observed in H a emission. Furthermore, 
the N-arm is well correlated with a northern dust lane char- 
acterized by a red V-I color, as derived from the HST images 
of the galaxy (Fig. 0). In contrast, the S-arm can be hardly 
identified in the V-I color image of Fig. [6] This overall north- 
south asymmetry in the color image can be partly explained 
if we assume that the northern side of the disk is the near 
side of the galaxy. However, the differences in the morphol- 
ogy of the N-arm and the S-arm, revealed by CO and H a , are 
to a large extent intrinsic and thus cannot be explained by a 
projection effect. 

- The outer ring: further out, north, there is a chain of emission 
clumps detected in CO(1-0) (hereafter referred to as outer 
ring or O-ring). This component, only detected in the 1-0 
line, has no clear southern equivalent. As shown in Fig.[6j the 
O-ring is well correlated with a chain of gas spurs character- 
ized by a red V-I color. Furthermore, the O-ring is also de- 
tected in H a emission (see Fig.[12]and discussion in Sect. [7]); 
this morphology, first described by Pogge ( 1989) as unusual, 
mimics a one-sided ellipsoidal ring. 

In summary, the bulk of the gas response in the disk of 
NGC 4579, as traced by the PdBI CO maps described above, 
follows the expected gas flow pattern induced by the bar poten- 
tial of the galaxy. In particular, the inner spiral arms appear as 
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Fig. 7. a)(Left panel) Same contours for the CO(2-l) line as in Fig. [6] overlaid on the OIII line color image from HST of Pogge et 
al. (12000) (color scale) observed in the nucleus of NGC 4579. Complexes N, N' , E' and E are highlighted. b)(Right panels) CO(2-l) 
line spectra observed towards complexes N, N', E' and E. Velocities are relative to v sys . 
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Fig. 9. a) Position-velocity diagram of CO(1-0) along the mayor axis of NGC 4579 (PA=95°). Contour levels go from -1.5<x, 1.5<x, 
3<x, 5<x, 8<x to 24<T in steps of 4<x where the 1-sigma rms <x=2.3 mJy beam -1 , b) Same as a) but for CO(2-l) line with levels going 
from -1.5<x, 1.5cr, 3<r, 5<r, Icr to 15<r in steps of 2<x where the 1-sigma rms cr=4.6 mJy beam -1 , c) Position- velocity diagram of 
CO(1-0) along the minor axis of NGC 4579 (PA=5°). Contour levels are -1 .5cr, 1 .5cr, 3cr and 5cr. d) Same as c) but for the CO(2-l) 
line. Contour levels are -1.5<r, 1.5<r and 3<r. In all panels velocities are relative to the systemic velocity (i^/=1470 km s" 1 ) and 
both Ax r and Ay r offsets along the major and minor axis, respectively, are relative to the AGN locus. The position of the N and E 
complexes and that of the O-ring are highlighted. 



two offset gas lanes which are leading with respect to the bar. gion in the nucleus, i.e., leaving room for two ILRs. In terms of 
This geometry suggests that there is spatially extended ILR re- bar orbit structure, this result can be interpreted as an indication 
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Fig. 8. a) CO(l-O) (Upper panel) and CO(2-l) (Lower 
panel) isovelocities contoured over false-color velocity maps. 
Velocities span the range (-225 km s" 1 , 225 km s" 1 ) in steps 
of 25 km s" 1 . Velocity scale is relative to v sys . Isovelocities have 
been derived using a 4-cr clipping on both data cubes. The AGN 
position is marked with a star. 

that there are stable stellar orbits whose major axis is perpen- 
dicular to the bar major axis (%2 orbits)(e.g., Athanassoula 1992; 
Buta & Combes [T9961 Regan & Teuben 12003]) . However, we 
have identified in the CO maps departures from the purely m-2 
point- symmetric gas flow pattern dictated by the bar. The re- 
ported deviations suggest that a superposed lopsided (m=l) in- 
stability may be at work. Firstly, there are indications of lopsid- 
edness in the inner CO spiral arms, echoed by other tracers of the 
ISM. The O-ring is an asymmetric pattern that could be the relic 
of a bar resonance, which is currently being depopulated. Closer 
to the AGN, the overall morphology of the central disk also sug- 
gests that the m-2 point- symmetry of the gas flow breaks up at 
r <200pc letting lopsidedness take over. 

5. CO line ratios 

The value of the 2-1/1-0 integrated intensity ratio (R21) is ob- 
tained after convolving the 2-1 map to the lower resolution of the 
1-0 map, both maps including short spacings. R21 ranges from 
-0.4 to -1 inside the observed region (in T m t km s" 1 units); 
the average ratio is -0.6. This ratio is close to the canonical 
value observed in spiral disks and is typical of optically thick 
molecular gas clouds (e.g., Braine & Combes 119921 Garcia- 
Burillo et al. 119931 ). There is a radial trend in the excitation 



of the gas, as measured by the line ratio: close to the AGN, at 
r <200 pc, R21 -0.8-1 (over the central disk), whereas R21 -0.4- 
0.6 at r>500 pc (over the S-arm, the N-arm and the O-ring). 
Hints of a similar trend are found in other NUGA targets (e.g., 
Garcfa-Burillo et al. 2003b) and can be interpreted as evidence 
for external heating of molecular clouds by X-rays in the vicinity 
of the AGN (Baker et al. 12003]) . 

6. Molecular gas kinematics 

6.1. Rotation curve and dynamical mass 

We show in Fig. [8] the isovelocity contour maps obtained from 
the CO(1-0) and CO(2-l) data in the nucleus of NGC 4579 (ap- 
plying a 4-cr clipping). The kinematics of the molecular gas are 
consistent with those of a spatially resolved rotating disk. The 
sense of the rotation of the gas in the galaxy plane is coun- 
terclockwise assuming the northern side is the near side (see 
Sect. g]). 

An estimate of the CO rotation curve (v rot ) can be obtained 
from position-velocity (p - v) diagrams taken along the ma- 
jor axis of NGC 4579 (Fig. [9]), assuming that circular motions 
dominate the gas kinematics. We have fitted the major axis po- 
sition angle as PA=95±10°; this value maximizes the line-of- 
sight velocity gradient within r=2 kpc. Furthermore, this de- 
termination is consistent within the errors with previous find- 
ings based on H a and HI kinematics (Guhathakurta et al. 1 1988t 
Chemin et al. 2006J Haan et al. I2008ab . Terminal velocities 
(Pterm) are derived by fitting multiple Gaussian profiles to the 
spectra across the CO major axis and selecting at each position 
the best estimate for v term , taking into account velocity disper- 
sion. During the fitting process we have masked the CO spectra 
between 0.5"(50 pc) and 2.5"(250 pc) east of the AGN because 
strong non-circular motions in this region make the gas appar- 
ently counter-rotate around the nucleus (see Fig. [9] and discus- 
sion in Sect. 16.21) . For the same reasons, the data points from the 
O-ring have not been included in the fit. The velocity centroids, 
corrected for inclination /=36°, and referred to v=v sys , give an 
estimate of v rot for each offset along the major axis. A simple 
rotation curve model can be fitted using the CO data from both 
sides of the major axis. We choose a rotation law characteristic 
of a logarithmic gravitational potential, described by: 



^roi\' / r— v v 

A/r 2 + r 2 

V ' ' max 

where v max is the maximum rotation velocity and r max is a 
radial scale whose value controls the steepness of v rot . To ex- 
tend v rot to the outer radii of the galaxy we have included in 
the global fit th e HI da ta points of Guhathakurta et al. (fT988t 
and Haan et al. (2008a). The best fit gives v max =295 km s" 1 and 
tw=530 pc (Fig. [10]). The v rot curve shows a gradual increase to 
r -0.65 kpc, where v rot -240 km s" 1 . From r -3 kpc to the outer 
disk traced by HI, up to r~ 12 kpc, v rot remains roughly flat at 
-295 km s" 1 . The mass inside a certain radius r can be inferred 
using M(r)=C x r x v 2 r JG, where G is the constant of gravity, 
M(r) is the mass inside a sphere of radius r, and C is a constant 
varying between 0.6 and 1, depending on the disk mass model 
assumed. If we take a value of C=0.8 intermediate between the 
values appropriate for spherical (1) and flat disk (0.6) distribu- 
tions, M(r -2 kpc)~4xl0 10 MQ. This implies that the molecular 
gas mass fraction inside r =2 kpc is -2%. This fraction decreases 
to <0.6% at r =50 pc. 

From the best fit of v roU we have derived the frequency 
curves, shown in Fig.[T0j Q, Q ± k/2 and Q - k/4. These curves 
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Fig. 10. a)(Upper panel) Rotation curve v rot and angular speed 
curve Q fitted from the CO and HI data of NGC 4579 (tri- 
angles are data points). A value for the bar pattern speed 
Qbar ~50 km s" 1 kpc -1 puts the corotation radius of the bar 
at rcR ~6 kpc. Similarly, a corresponding pattern speed value 
of Qoval -270 km s" 1 kpc -1 puts the corotation radius of the 
oval at ~1 kpc. b)(Lower panel) Frequency curves £1 ± k/2 
and Q - k/4. If Qbar ~50 km s" 1 kpc -1 , the Ultra Harmonic 
Resonance (UHR) of the bar lies at r UH R -3.8 kpc. We estimate 
the existence of an extended ILR region with two ILRs lying 
at ruLR -500 pc (inner ILR=iILR) and r oILR ~1.3 kpc (outer 
ILR=oILR), respectively. With Qoval ~270 km s" 1 kpc -1 the 
UHR of the oval is at -200 pc. 



can be used to estimate the position of the main resonances in the 
disk in the context of the epicyclic approximation as discussed 
in Sect. [71 where we compare the morphology of the molecular 
disk to the distribution shown by other tracers of the ISM and 
the stellar structure in the disk of NGC 4579. 



6.2. Streaming motions 

The kinematics of molecular gas in the central r~2 kpc of 
NGC 4579 reveal the presence of non-circular motions that 
mod ulate the overall rotation pattern of the disk described in 
Sect. 16.11 Deviations are detected with uneven strength over the 
inner spiral arms (N-arm and S-arm), the outer arm (O-ring), 
and, most noticeably, over the central disk, as described below: 

- The central disk: the radial velocities of the gas in the cen- 
tral disk deviate by up to >100kms _1 from the expected pat- 
tern of circular rotation. CO emission from the E complex, 
mostly prominent in the 2-1 line, is detected at negative ra- 
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Fig. 11. a) We overlay the HI intensity map obtained with 
the VLA (in contours from 0.04 to 0.28 in steps of 
0.03 Jy km s" 1 beam -1 ) on the ^-band WHT image (color 
scale) obtained for the nucleus of NGC 4579; both images have 
been deprojected onto the galaxy plane. Units on X/Y axes 
(AXg/AYg) correspond to arcsec offsets along the major/minor 
axes with respect to the AGN. The orientation of the large-scale 
bar-BAR-is shown, b) Same as a) but we overlaid the CO(l- 
0) PdBI intensity map (contour levels=0.2, 0.5, 1, 1.7, 2.6, 3.7 
and 5 Jy km s" 1 beam -1 ), c) Same as a) but for the CO(2-l) 
PdBI intensity map (contour levels=0.3, 0.6, 1 to 9 in steps of 
1 Jy km s" 1 beam -1 ). 



dial velocities ( — 70 kms -1 ) that axe forbidden by circular 
rotation at this location of the disk (see Figs.[7]and [9]). If me 
gas flows inside the plane, this would imply that the gas at E 
is counter-rotating at a deprojected speed of f;~120kms _1 . 
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Complex N, which lies close to the kinematic minor axis 
of the galaxy, also shows non-circular motions: instead of a 
zero average velocity, we measured — 50 kms -1 (see Figs. [7] 
and [9]); this can be interpreted as gas moving outward within 
the plane. The CO kinematics observed in complexes N' 
and E' indicate instead positive radial velocities for the gas: 
~+75 kms -1 at N' and ~+25 kms -1 at E'. Altogether, the ra- 
dial velocities measured along a line intersecting complexes 
N, N', E' and E show a gradient that is reversed with respect 
to that expected if circular rotation prevailed. 

- The inner spiral arms: the isovelocity maps of Fig. [8] 
reveal a steep transverse velocity gradient across the 
S-arm. An inspection of the major axis p-v dia- 
grams of Fig. [9J2,/? indicates that the gas is strongly 
slowed down to zero rotation velocities when it en- 
counters the S- arm at Ax'~+3". Downstream the S- 
arm (at Ax'~+7"), gas is accelerated: rotation velocities 
v ro t=(v - ^)/sin(/)~170/sin(36)~290 km s" 1 are reached 
when the gas leaves the arm. A similar but more shallow 
transverse gradient is detected across the N-arm. 

- The outer ring: Figure [8] shows that the O-ring is also asso- 
ciated with peculiar motions. The radial velocities measured 
in adjacent regions of the N-arm and O-ring can differ by up 
to 150 km s" 1 . This is illustrated in Fig. [9] which shows the 
p-v plots along the major and minor axis of the galaxy. 

The streaming motions detected over the inner spiral arms 
are in agreement with the expected kinematic signatures of the 
stellar bar on the gas flow within the ILR region, and thus, well 
inside corotation of the bar (e.g., see Canzian 1993 ). However, 
the kinematic differences between the N-arm and the S-arm 
point to the influence of a superposed lopsided instability that 
could also be at play in the central disk component. The kine- 
matic signature of lopsidedness is also found in the O-ring. Non- 
circular motions let the O-ring appear as a separate component 
decoupled from the inner gas response (N-arm, S-arm). 

6.3. An outflow of molecular gas in NGC 4579? 

Alternatively, the overall morphology and velocity pattern of 
the central disk could be qualitatively explained by a scenario 
where the gas is flowing away from the nucleus in extra-planar 
motions, as mentioned by GB05. The OIII line and UV contin- 
uum maps obtained with HST can be used to tentatively iden- 
tify a one-sided fan-shaped morphology in the emission com- 
ing fro m the c entral r ~ 100 pc of NGC 4579 (Pogge et al. 120001 
Contini 2004 ), suggestive of an AGN-powered outflow of hot 
gas. The morphology of the neutral gas disk at r<150 pc, evi- 
denced by the V-I color image of the galaxy, would indicate that 
the gas is piling up in a shell. The shell could be a signature of 
gas entrainment by the outflow (see Fig. [7]). The expanding shell 
scenario can explain the kinematics of the N and E complexes as 
a signature of out-of-the-plane motions projected along the line 
of sight towards the observer. The positive radial velocities mea- 
sured towards complexes N' and E' can be explained if out of 
the plane motions are projected away from the observer. In this 
scenario AGN feeding would be presently quenched on these 
scales. 

The outflow scenario for NGC 4579 is not flawless, however. 
Firstly, based on the current data, the putative outflow of hot 
gas in NGC 4579 does not have the mass and energies typically 
required to produce a corresponding outflow in the neutral gas 
(e.g., see the case of NGC 4528 discussed by Krause et al. l2007b . 
Furthermore, in NGC 4579 there is no evidence of a radio jet, 
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Fig. 12. Same contours as in Fig. [6] overlaid on the H a image 
from HST (color scale) observed in the nucleus of NGC 4579. 

which should be the driving agent of the outflow. Finally, there 
is no information on the kinematics of the ionized gas that may 
be indicative of outflow motions. On the whole, the outflow sce- 
nario can be considered as unlikely and merely speculative. 

7. Comparison with other tracers 

In this section we compare the CO maps of NGC 4579 with the 
distribution shown by other tracers of the ISM and with the stel- 
lar structure of the galaxy. This information is used to estimate 
the location of the resonances in NGC 4579, a key to the inter- 
pretation of the gas-flow in the disk of the galaxy. 

7.1. The large-scale bar 

T he 2.2/i m ^T-band image of NGC 4579 published by Knapen et 
al. l2003l shows a prominent stellar bar in the disk. The bar is ori- 
ented along PAb ar ~56±5° (measured from north in the plane of 
the sky). Once deprojected onto the plane of the galaxy (assum- 
ing PA=95° and /=36°), the position angle of the bar (PA' bar ) 
measured from the major axis becomes: tan( J PA^ ar )=tan( J PA- 
PAb ar )/cos(i), which implies PA' bar -45°, as shown in Fig. [TT] 
An isophotal analysis made on the K-band image of the galaxy 
gives an approximate diameter- size for the bar of D -10-12 kpc. 
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Fig. 13. a)(Upper panel) Superposition of the HI intensity map 
(black contours), the 3.6 jim Spitzer/IRAC image in contour lev- 
els 1.1, 1.3, 1.6, 2, 3, 4, 6 and 11 MJy/sr (white contours) and the 
8 yum Spitzer/IRAC image (grey scale). All images have been de- 
projected onto the galaxy plane. b)(Lower panel) Same as a) but 
HI intensity map (in contours) superposed on the FUV GALEX 
image (grey scale). 
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Fig. 14. a)(Upper panel) The 12 CO(2-l) integrated intensity 
map obtained with the PdBI (same contours as in Fig. [6]) is over- 
laid on the UV image of the central region of NGC 4579 obtained 
with the ACS/HRC in the F330W band (in grey scale). b)(Lower 
panel) The ^T-band contours are overlaid on the UV-image of the 
galaxy. We highlight the orientation of the bar and the oval. All 
images have been deprojected onto the galaxy plane. 



This coincides with the value derived from a similar analysis 
made on the 3.6/mi Spitzer/IRAC image (Fig. [13]). This bar size 
is compatible within the errors with the value calculated from the 
Fourier decomposition of the potential, which favors a slightly 
larger bar (D >14 kpc), as discussed in Sect. 18.2. II 

In the following, we will suppose that the corotation res- 
onance lies at the end of the bar: vcr ~6±1 kpc, i.e., al- 
lowing for a 20% uncertainty in vcr . This is a plausible 
scenario for the well-developed although relatively weak bar 
of NGC 4579 (see Sect. [8]). An estimate of the bar pattern 
speed (Qbar) is thus obtained using the D curve derived from 
CO if we rely on the epicyclic approximation; this yields 
&bar -50+10 km s _1 kpc _1 (See Fig. [TO]). The two ILRs lie ap- 
proximately at r -500 pc (inner ILR=iILR) and r ~ 1 .3 kpc (outer 
ILR=oILR) for £l BAR -50 km s^kpc" 1 (Fig. [TO]). This value of 
Qbar allows the bar to have an extended ILR region, a require- 



ment imposed by the observed distribution and kinematics of 
molecular gas in the central r ~2 kpc disk. In the context of this 
approximation we can put an upper limit to the pattern speed 
of the bar of NGC 4579: £l BAR ~60 km s^kpc" 1 ; this value of 
Qbar would impose the existence of a single ILR at -800 pc, a 
scenario incompatible with observations, which require the ex- 
istence of two ILRs. With Q#a#~50 km s _1 kpc _1 , the O-ring 
would correspond to the oILR, a resonance which is currently 
being depopulated due to gas inflow down to a region inside the 
ilLR of the bar (see discussion in Sect. 18.21) . The latter explains 
the molecular gas concentration in the inner spiral arms (N-arm 
and S-arm), formed as a mixture of X2 and xi precessing or- 
bits. On its way to the nucleus, molecular gas is forming stars in 
the inner spiral arms, as evidenced by the detected H a emission, 
more intense in the case of the N-arm (see Fig. [12]), as discussed 
in Sect. EH 
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7.2. The nuclear oval 

The ^-band image of NGC 4579 shows evidence of a previously 
unnoticed oval distortion present in the inner r -200 pc (Tig.[T4b. 
This nested structure was not detected in the /-band HST image 
originally used in the first gravity torque analysis of NGC 4579 
(GB05). The oval is oriented along PA' oval -70° (measured from 
+X axis in the plane of the galaxy), and therefore is not aligned 
with the large-scale bar. Misalignment can be taken as evidence 
of decoupling of the nested structure. The reality of this distor- 
tion is confirmed by the analysis of the Fourier decomposition of 
the potential described in Sect. [8] Fig. [T4l shows a zoomed view 
of the central r -200 pc-region. The oval appears encircled by 
an off-centered ellipsoidal ring detected in the HST UV image 
of the galaxy, first discussed by Contini (2004). The UV ring 
has an average radius r -150-200 pc and has been recently de- 
scribed by Comeron et al. (2008) as an ultra-compact nuclear 
ring. This ring can be interpreted as a resonance of the oval dis- 
tortion. The best guess assigns the Ultra Harmonic Resonance 
(UHR) of the oval to the UV ring: if the pattern speed of the oval 
is tlovAL -270+10 km s _1 kpc _1 , the UHR of the oval should 
lie at -200 pc and its corotation at -1 kpc (Fig.[Kj). With this 
choice we can explain the decoupling of the oval with respect to 
the large-scale bar, as the corotation of the oval would be close 
to the oILR of the large-scale bar. This resonance overlap is a 
prerequisite to the decoupling of nuclear bars/ovals (e.g., Garcia- 
Burillo et al. [T9981 Hunt et al. [2008b . 

As discussed in Sect. 18.21 the oval perturbation is key to fa- 
voring gas fueling to the AGN at present. More particularly, the 
decoupling scenario explains the depopulation of the oILR of 
the large-scale bar, under the combined action of the bar and the 
nuclear oval. A similar time sequence is nicely illustrated by the 
numerical simulations adapted to follow the nuclear bar decou- 
pling in the double-barred galaxy NGC 2782, discussed by Hunt 
et al. (120081) . 

7.3. The gas distribution in the outer disk (r>4 kpc) 

The 21cm HI map of NGC 4579 gives a sharp view of the gas 
distribution in the outer disk of the galaxy (from r -4 to 15 kpc) 
(See Figs. and [13]). The 8/mi Spitzer/IRAC and the FUV 
GALEX images provide a detailed picture of the star formation 
in the disk. For a Sb galaxy like NGC 4579, >70% of the total 
emission at 8/mi is of non-stellar origin and is produced by warm 
dust heated by star formation activity (e.g., Pahre et al. 120041) . 
The HI distribution in the disk shows a central depression in the 
inner r -4 kpc, mimicking a ring. The HI pseudo-ring is partly 
incomplete and has an average diameter of -12 kpc. At close 
sight, the HI pseudo-ring appears to be the superposition of two 
winding spiral arms. The inner segments of the spiral arms are 
better identified in the 8/zm Spitzer/IRAC and the GALEX im- 
ages; the 8//m and FUV arms neatly delineate the inner edge of 
the HI pseudo-ring. Further out, the arms unfold in the disk over 
180° in azimuth and end up forming the pseudo-ring detected in 
HI. Most remarkably, the spiral arms are not connected to the 
two bar ends located at (r,<D)=(6 kpc, 45°) and (6 kpc, 225°). 
Instead, the two-arm spiral structure starts at a different azimuth 
(-90° from the bar ends) and, also, at a significantly smaller ra- 
dius in the disk: (r,0)=(4 kpc, 135°) and (4 kpc, 315°). 

On the whole, the particular geometry described above sug- 
gests that the bar and the spiral are decoupled and thus do not 
share a common pattern speed. It is accepted that spirals and 
bars can rotate at differing angular speeds in barred spiral galax- 
ies, a scenario first advocated by Sellwood & Sparke ( 198S]) who 



contemplate that, under decoupling, outer spirals can rotate at 
lower £l p . As recently shown by numerical simulations (e.g., 
Combes [2007 1 and references therein), decoupling is more likely 
in the case of weak bars, of which the NGC 4579 's bar is a good 
example. As discussed in Sect. 18.21 this situation favors gas fu- 
eling to the central regions, as the corotation barrier imposed by 
the bar could be thus overcome. The morphology of the outer 
disk revealed by HI in NGC 4579 suggests that the neutral gas is 
currently piling up in an inner ring located well inside the coro- 
tation of the bar. This scenario is confirmed by the coincidence 
of the inner edge of the HI pseudo-ring with the predicted posi- 
tion of the UHR of the bar, around r UHR -3.8 kpc, as illustrated 
in Fig. [TO] 

7.4. The star formation rate in NGC 4579 

The star formation pattern in the inner disk (r <2 kpc) 
of NGC 4579 have been previously described by Pogge et 
al. (12000b and Contini (120041) . based on HST optical and UVim- 
ages (see also descriptions by GB05 and Comeron et al. 12008b . 
Most of the H a emission in the region 200 pc< r <2 kpc comes 
from a two-arm spiral structure that shows a geometry similar 
to that of the CO inner spiral arms, described in SectH (see 
Fig. [12]). In addition, part of the H a emission is related to the 
CO outer ring. Closer to the AGN (r <200 pc), the bulk of 
the H a emission is seen to come from a slightly off-centered 
disk. At close sight, several subcomponents can be identified in 
this disk. There is a point-like nuclear source of ~0.27"-size, 
located at the dynamical center of the galaxy; the emission of 
this component has been in terpret ed as due to photo-ionization 
by the AGN (Maoz et al. 2005). Furthermore, the OIII line 
and UV maps reveal also a one-sided fan-shaped morphology 
in the emission in the inner r -100 pc of NGC 4579 (Pogge et 
al. [20001 Contini [^0i). On larger scales, the HST UV image 
of the galaxy, shows an off-centered ring at r -150-200 pc (see 
Sect. 17.21 and Comeron et al. 2008). As described in Sect. 17.31 
star formation in the outer disk, as traced by the UV continuum 
emission seen by GALEX, arises from a pseudo-ring character- 
ized by a geometry similar to that of the HI pseudo-ring. 

In order to quantify the amount of gas transformed into 
stars, we have used the HST H a and GALEX UV images de- 
scribed above to calculate the corresponding star formation rates 
(SFR) on different regions throughout the disk of the galaxy. 
SFR values are derived using the empirical formulae given by 
Kennicutt (1998). We first estimate the SFR in the inner disk; 
more precisely, we define S¥R inner as the integrated SFR inside 
the region 300 pc< r <2 kpc. To correct for extinction, and thus 
have a more accurate evaluation of S¥R inner , we have used the 
CO(1-0) PdBI map to derive the H2 column density map, using 
the CO-to-H2 conversion factor of Solomon & Barrett (11991b : 
X=N(H 2 )/Ico=23xl0 20 cm" 2 K" 1 km" 1 s. The A y map is then 
derived using the commonly-used N(H2)-to-A y conversion fac- 
tor of Bohlin et al (j 19781) : N(H 2 )/Ay=10 21 cm- 2 mag _1 . Adopting 
a screen-geometry for the obscuring dust in front of the HII re- 
gions, we must scale the A v map by a factor 1/2, if HII regions 
are assumed to lie in the mid-plane of the galaxy. This scaled ex- 
tinction map, with A v typically ranging from -3 to -15, is used 
to derive A Ha , following Cardelli et al. ([T989b : A# a =0.828xA y . 
We derive SFR inner -0.10 Mq yr" 1 , confirming that the nuclear 
starburst of NGC 4579 is only moderate. In a second step, we use 
S¥R inner to calculate the SFR throughout the disk of the galaxy 
by scaling the integrated fluxes measured in the GALEX UV im- 
age at the corresponding radii. In particular, we define S¥R 0Uter 
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Fig. 15. Strength (Q\, Q2 and gr) and phases (<p\ and $2) of m = 1 and m = 2 Fourier components of the stellar potential derived 
from the K-band WHT image inside r~15 kpc (left panels: sl, b) and inside r~2 kpc (right panels: c, d). We highlight the radii 
where the OVAL and the BAR perturbations prevail at r <2 kpc. In addition, Q3, and Q4 values are plotted in panel a to illustrate 
the predominance of the bar mode in the disk up to r ~7-8 kpc. 



as the SFR integrated inside the region 4 kpc<r<10 kpc. We 
derive SFR OM?er ~1.30 Mq yr -1 . 

As such, SFRinner and S¥R 0Uter contain the bulk of the total 
SFR of NGC4579 estimated as S¥R tota i -1.70 M Q yr" 1 . These 
SFR are compared to the corresponding inflow rates driven by 
gravity torques in these regions, relevant to understand the over- 
all gravity torque budget, as discussed in Sect. 18.2.31 



8. Gravity torques 

8.1. Methodology 

Firstly we have derived a good representation of the stellar po- 
tential for the disk of NGC4579. To do so we have used the 
^-band image of the galaxy published by Knapen et al. (120031) , 
assuming that this image is virtually free from dust extinction 
or stellar population biases. This working hypothesis better ap- 
plies to the ^T-band image than to the /-band image originally 
used by GB05 to derive the gravity torque budget in NGC 4579. 
The ^T-band image has been deprojected according to the an- 
gles PA=95 and i=36 °, assumed in Sect.[T] Figure. [TT1 shows 
the overlay of the deprojected gas distribution as derived from 
the CO and HI maps with the deprojected K-band image of 
NGC 4579. 

To reproduce the vertical mass distribution, the K-band im- 
age was convolved using an isothermal plane model of constant 
scale height, equal to ~ 1/1 2th of the radial scale-length of the 
image. We also adopted a constant mass-to-light (M/L) ratio, 
obtained by fitting the rotation curve to the CO and HI data. 
The potential-®(7?, #)-was then calculated by a Fourier trans- 
form method and decomposed in the different m-modes as fol- 
lows: 



®(R, 6) = O (R) + V <b m (R) cos(m6 - <j> m (R)) 



(2) 



where $> m (R) and <f> m (R) are the amplitude and phase of the m- 
mode, respectively. 



According to Combes & Sanders (1981), the strength of the 
m-Fourier component is defined as: Q m (R) as 



Q m (R) = m^ m /R\F (R)\ 



(3) 



Similarly, the strength of the total non-axisymmetric pertur- 
bation is obtained by: 



Qt(R) 



F™ X (R) k{ 



l(dtyR,0)\ 



36 



Fo(R) 



dR 



(4) 



where F™ ax (R) represents the maximum amplitude of the tangen- 
tial force over all 6 and Fo(R) is the mean axisymmetric radial 
force. 

Figure[l5] shows the description of the gravitational potential 
of NGC 4579 defined by the Q/=i,2, Qr, and 0/=i,2 curves, used 
in the discussion of Sect. 18.2.11 

The forces per unit mass (F x and F y ) are calculated from the 
derivatives of 0(7?, 6) at each pixel; the gravity torques per unit 
mass-£(x, z/)-can be then computed by: 



t(x, y) = xFy-y F x 



(5) 



The sense of the circulation of the gas in the galaxy plane 
(counterclockwise in NGC 4579; see Sect. [6]) defines whether 
the sign of t(x, y) is positive (negative) if the torque accelerates 
(decelerates) locally the gas in the disk. Note that the specific 
torque field t(x, y) does not depend on the present distribution 
of the gas. However, we use the t(x, y) map to infer the angular 
momentum variations by assuming that the gas column density 
(N(x, y)) derived from the CO line maps (in the inner disk) and 
from the HI map (in the outer disk) is a good estimate of the 
probability of finding gas particles at this location at present. The 
torque field is then weighted by N(x, y) at each pixel to derive the 
time derivative of the angular momentum surface density in the 
galaxy plane, dL s (x, y)/dt=N(x, y) x t(x, y). Our approach is sta- 
tistical since we do not follow particles along individual orbits, at 
best estimated through model fitting (e.g., see Boone et al. 120071) , 
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Fig. 16. a)(Upper panel) We overlay the HI contours with the the 
map of the (dimensionless) effective angular momentum varia- 
tion, as defined in text. The derived torques change sign as ex- 
pected, if the butterfly diagram, defined by the orientation of 
quadrants I-to-IV, can be attributed to the action of the large- 
scale bar. b)(Middle panel) The same as a) but for CO(l-O). 
c)(Lower panel) The same as a) but for CO(2-l). Maps are de- 
projected onto the galaxy plane. 



in order to derive angular momentum variations. Instead, we av- 



erage over all orbits of gas particles at each position and account 
for the time spent by the gas clouds along all the possible orbit 
paths. Figure [16] shows the dimensionless version of the torque 
maps normalized by [| N(x, y) x t(x, y) \] max , as derived from CO 
and HI. Sect . 1 8 . 2 . 2l pro vides a detailed description of these maps. 

To estimate the net radial gas flow globally induced by the 
torque field in the galaxy disk, we have computed the torque per 
unit mass averaged over the azimuth 6 as a function of radius R, 
using N(x, y) as weighting function: 



t(R) 



f g N(x,y)x(x F y -yF x ) 



(6) 



By definition, t(R) is the azimuthally-averaged time deriva- 
tive of the specific angular momentum-L-of the gas, i.e., 
t(R)=dL/dt\o. Similarly to the torque maps, the sign of t(R), ei- 
ther positive or negative, defines whether the gas gains or loses 
angular momentum, respectively. Specifically, we estimate the 
AGN feeding efficiency by deriving the average fraction of the 
gas specific angular momentum transferred in one rotation (T rot ) 
by the stellar potential, as a function of radius, i.e., by the di- 
mensionless function AL/L defined as: 



AL 



dL 
dt 



1 



XT rt 



KR) 



xT r( 



(7) 



where L e is assumed to be well represented by its axisymmet- 
ric average, L e = R x v rot . The absolute value of L/AL deter- 
mines how much time (in T rot units) the stellar potential will 
need to transfer the entirety of the gas angular momentum. A 
small value of AL/L implies that the stellar potential is ineffi- 
cient at present. The AL/L curves derived from the CO and HI 
maps of NGC 4579, shown in Fig.[17j are used in the discussion 
of Sect. [8221 

Following the approach of GB05, we can estimate the gas 
mass inflow (- sign)/outflow (+ sign) rate driven by the stel- 
lar potential per unit length as a function of radius (in units of 
M yr _1 pc _1 in Figs.[T9b tolSTh) according to: 



d 2 M 
dRdt 



dL 
dt 
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e L 



x 2nR x N(x, y) \ e 



(8) 



where N(x, y) \ e is the radial profile of N(x, y) averaged over the 
azimuth for a radial binning AR, derived from HI, for the outer 
disk, and from CO, for the inner disk. 

Similarly, the inflow/outflow rates integrated out to a certain 
radius R can be derived as: 



dM _ y d 2 M 



xAR 



(9) 



Figs. lT9b to [2Tb display these integrated rates in units of 
M© yr" 1 . The mass transfer rates in NGC 4579 are discussed in 
Sect. 18X31 

8.2. The gravity torque budget in NGC 4579 
8.2.1 . The stellar potential 

To derive the stellar potential 0(7?, 6) used in the gravity torque 
calculation of NGC 4579, we have followed an adaptive ap- 
proach. ®(R,6) in the inner r~2 kpc disk has been derived 
using a grid of 256x256 pixels with 0.15" pixel -1 on the K- 
band image. This representation of 0(/?, 6) is the best adapted 
to derive the torque field on CO. On the other hand, a spatially 
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smoothed version of the NIR image, with 256x256 pixels and 
1.15" pixel -1 , has been used to calculate $>(R,6) in the outer 
disk up to r -15 kpc. The ensuing torque field on HI is derived 
from this representation of the stellar potential. 

Figure Q3] shows the Q values derived as a function of radius. 
The large-scale bar (m = 2) is clearly detected in the smoothed 
version of 0(7?, 6) (Figs.fTSb,^). The area of influence of the bar, 
determined by the dominance of the m = 2 mode over other sec- 
ondary modes extends up to r -7-8 kpc in the disk. For r >8 kpc, 
G1+G3+G4 ^62, an indication that the bar mode is not domi- 
nant on these scales of the outer disk. This puts an upper limit 
to the size of the bar (r<7-8kpc). As discussed in Sect. 17.11 
the isophotal analysis of the NIR images of the galaxy gives a 
smaller size for the bar: r -5-6 kpc. The finally adopted value 
for vcr =6+1 kpc is a reasonable compromise between the two 
determinations above and it accounts for the overall gas response 
in the disk. 

On the whole, the bar is moderately strong: the Q2 strength 
of the bar is maximum at r~3 kpc where it reaches a value 
of -0.1. Moreover, the phase of the bar 02 stays noticeably 
constant at a value — 0.8 radians, i.e., at an angle -45°, mea- 
sured counterclockwise from +X-axis in Figs.[TT] This also con- 
firms the value of PA' b derived from the isophotal analysis of 
Sect. 17.11 The inner representation of 0(7?, 6) also reveals the 
bar at r<2 kpc (Figs. [T5l:,<fh with a similarly constant phase 
and, not surprisingly, with an average lower strength, Q2 -0.02- 
0.03 (at these radii we are inside the peak of Q2). Of particular 
note, the high-resolution version of <D(7?, 6) unveils the oval dis- 
tortion with a size and an orientation similar to those derived 
in Sect. 17.21 The presence of a misaligned oval is betrayed by a 
significant -25° change in the phase of the m=2 mode. 02 gradu- 
ally shifts from the large-scale bar orientation, — 0.8 radians for 
200 pc< r <7 kpc (-45°), to —1.2 radians (-70°) at r -100 pc. 

There are discrepancies between the stellar potential ob- 
tained from the AT-band image of NGC 4579, and the represen- 
tation of <D(#, 0) obtained by GB05 from the /-band HST image 
of the galaxy. Firstly, the strong m = 1 mode detected in the 
/-band image is much weaker in the new evaluation of 0(7?, 0). 
As illustrated in Fig. 7 of GB05, the m = 1 mode seemed to 
overtake the potential strength of the bar at r<500 pc: Qi was 
seen to vary between -0.02 and -0.10, while Q2 -0.02. This is 
in contrast to the new result shown in Fig.[T5b: Q2 > Qi, except 
for the inner r -50 pc, i.e., below the limit of the spatial resolu- 
tion of the imageQ. This has consequences on the interpretation 
of the lopsidedness identified in the gas distribution and kine- 
matics probed by CO. Lopsidedness cannot be clearly attributed 
to the action of a true m = 1 mode in the stellar potential. As 
a second important difference, the oval identified in the K-bznd 
image was not present in the /-band image of GB05. The re- 
ported discrepancies in 0(7?, 6) can be attributed to the residual 
effects of extinction which are noticeable in the central r -2 kpc 
of the /-band image. 

Taken together, these differences are relevant to explaining 
the new version of the gravity torque budget in NGC 4579, de- 
scribed in Sect. 18.2.21 

8.2.2. Gas-flow time-scales in NGC 4579 

Figure [16] shows the 2D pattern of gravitational torques obtained 
throughout the disk of NGC 4579, from r - 15 kpc down to the 



4 The peak of Q\ at r~50-100 pc identified in Fig.rBb is an artifact 
created by the spatial smoothing performed to derive the large-scale 
representation of ®(R, 6), and as such is unreliable. 
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Fig. 17. The average fraction of the angular momentum trans- 
ferred from/to the gas in one rotation-JL/L-are plotted as a 
function of radius, as derived from the HI (upper panel), CO(l- 

0) (middle panel) and CO(2-l) (lower panel) maps of the disk 
of NGC 4579. The locations of BAR resonances (ilLR, oILR, 
CRbar, UURbar) and those of OVAL resonances (CRoval, 
UHRoval) are highlighted. We have blocked the central r -1 kpc 
in the dL/L plot obtained from the HI map, due to the spatial 
resolution limit of the HI data. For similar reasons, the central 
r - 100 pc and r -50 pc around the AGN are blocked in the dL/L 
plots obtained from the CO 1-0 and 2-1 maps, respectively. In 
the range 1.8 kpc> r >1.2 kpc, estimates of dL/L based on the 
CO(1-0) map are more reliable than those derived from CO(2- 

1) data (in grey), as the latter come from a region outside the 
PdBI primary beam at this frequency. The same display conven- 
tions are applied to Figs.[T9ltol2Tl 



inner r -50 pc around the AGN. The derived torques change sign 
following a characteristic 2D butterfly pattern. Quadrants I-to-IV 
will define hereafter the regions where the signs of the torques 
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as Fig. [T71 but here showing the results obtained 
CO maps with and without short spacing correc- 



cates a net gas inflow from the oILR down to a region inside the 
ilLR of the bar. As for the outer disk, the implied time- scales 
for angular momentum transfer are also -5-10 rotation periods. 
This result is in agreement with our first calculation published by 
GB05 (see Fig. 1 1 of GB05) regarding the sign of the torques. 
However the gravity torques shown in Fig [TTb , c are compara- 
tively stronger at the same radii, in particular those contributing 
to gas inflow. The changing picture can be attributed to the new 
0(7?, 6) used in this work but, also, to the new N(x, //)-values 
that now include the short spacing correction for CO. This is il- 
lustrated in Fig. [18] where we compare the dL/L curves obtained 
with the new <D(7?, 6) used in this work, but with two different 
versions of N(x, y), derived from the PdBI-only CO data (used 
by GB05), and from the short spacing corrected PdBI+30m data 
(this work), respectively. In summary, we find that the gas in 
the inner disk is funneled from the oILR to a region inside the 
ilLR at r -300 pc. There is a barrier of positive torques between 
r~200pc (i.e., close to the UHR of the oval perturbation) and 
r~300pc(Fig[T7k c). 

Closer to the AGN (r <200 pc), the bulk of the CO emission 
is concentrated north of the nucleus (complex N) in the trailing 
quadrant of the large-scale bar I(+), where stellar torques are 
positive. However, the rest of the central disk complexes from 
r -50 pc to r ~150pc, identified in the CO(2-l) map, feel nega- 
tive torques due to the combined action of the large-scale bar and 
the inner oval. The two m = 2 modes act in concert to produce 
net gas inflow down to the spatial resolution of our observations 
(see FigfTTb). This result is different to that reported in GB05, 
where the evidence of AGN fueling at r <200 pc was lacking. 
The implied time- scales for angular momentum transfer are ex- 
tremely short: ~l-3 rotation periods. 



driven by the dominant perturbation of the stellar potential are 
expected to be constant. 

Assuming that the gas is rotating counterclockwise, Fig.fT6b 
shows that, although HI emission is distributed along the four 
quadrants, the majority of the gas is seen to lie along the lead- 
ing edges of the large-scale bar where torques are negative 
(quadrants II(-) and IV(-)), especially on intermediate scales 
(4 kpc<r<8 kpc). Figure [TTh shows that the azimuthally av- 
eraged torques are systematically negative from r -10 kpc down 
to r ~4 kpc. This result suggests that the corotation barrier of the 
bar in NGC 4579 (at ~6 kpc) has been overcome due to secular 
evolution processes. If the gas response was the canonical spiral 
pattern coupled to the bar, the corotation would separate in the 
disk the regions where torques are positive (from corotation to 
the Outer Lindblad Resonance (OLR)) from those where torques 
are negative (from corotation to the oILR). Positive torques 
would tend to form a ring at the OLR. However, a decoupling of 
the spiral allows the gas to efficiently populate the UHR region 
inside corotation and thus produce net gas inflow on intermedi- 
ate scales in the disk. A similar scenario has been suggested to 
explain the gravity torque maps of some o f the NUGA galax- 
ies recently analyzed by Haan et al. (2008b). The implied time- 
scales for angular momentum transfer in NGC 4579 are -5-10 
rotation periods. On its way to the nucleus, the gas stops in the 
outer disk close to the UHR of the bar, where torques become 
negligible (FigfTTb). 

In the inner disk, Fig. [T6b , c shows that the bulk of the 
CO emission along the inner spiral arms lies where torques 
are strong and negative (quadrants II(-) and IV(-)). The corre- 
sponding azimuthally averaged torques, shown in FigfTTb, c, are 
mostly negative from r~1.3 kpc down to r~300 pc. This indi- 
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Fig. 19. a)(Upper panel) We represent the radial variation of the 
mass inflow(-) or outflow(+) rate of gas per unit radial length 
(d 2 M/dRdt) in the nucleus of NGC 4579 due to the action of 
stellar gravitational torques on the gas distribution as derived 
from HI data. Units are Mq yr _1 pc _1 . b)(Lower panel) We plot 
the mass inflow/outflow rate integrated inside a certain radius r 
(dM/dt) in M© yr" 1 . 
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Fig. 20. Same as Fig[l9]but here based on the gas distribution as 
derived from CO(1-0) data. 
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Fig. 21. Same as Fig[l9]but here based on the gas distribution as 
derived from CO(2-l) data. 



8.2.3. Mass transfer rates in NGC 4579 

In terms of mass inflow/outflow rates, the feeding budget inte- 
grated out to a certain radius is clearly negative in the outer disk 
of NGC 4579 (r >4 kpc). As shown in Fig.[l9l dM/dt goes from 
virtually zero inflow at r ~4 kpc to — 0.3 Mq yr" 1 at r ~10 kpc. 



The minimum of d 2 M/dRdt, indicative of the largest inflow rate 
per unit length, is reached at r~5.5 kpc, a position close to the 
corotation of the bar. As argued in Sect. 18.2.21 this paradoxical 
result can be explained if the gas response in the outer disk is 
decoupled from the bar potential. The inflow rate integrated on 
these scales (dM/dt[4 kpc< r <10 kpc] — 0.3 M© yr" 1 ) is a fac- 
tor ~4 lower than the gas feeding needs implied by the current 
SFR integrated on this region (SFR 0M?er ~1.3M© yr -1 , as derived 
in Sect. 17.4b . Note, however, that on these scales, the dM/dt esti- 
mated from HI is a lower limit to the total inflow rate, as it does 
not include the contribution from molecular gas at these radii. 

In the inner disk, the integrated mass feeding budget is sys- 
tematically negative from r ~2 kpc down to r -400 pc. As shown 
in Figs. [301 and ED dM/dt goes from — 0.5 M© yr -1 at the 
ilLR to — 1.5 M0 yr -1 at the oILR. This is a factor ~5 larger 
than the corresponding dM/dt inflow rates measured for the 
outer disk (r >4 kpc). The minimum of d 2 M/dRdt, indica- 
tive of the largest inflow rate per unit length, is reached in the 
vicinity of the ilLR. The inflow rate integrated on these scales 
(dM/dt[400 pc< r <2 kpc] —1.5 M Q yr" 1 ) is an order of mag- 
nitude higher than the gas feeding rate required to sustain the 
current SFR in this region (S¥R inner -0.10 Mq yr -1 , as derived 
in Sect. 17.4b . The time-scale for gas consumption by star for- 
mation is significantly larger than the dynamical time- scale for 
inflow in the inner disk at present. 

Closer to the AGN (r <200 pc), the integrated mass budget 
derived from the CO(2-l) map shows negative torques down to 
r ~50 pc. The minimum of d 2 M/dRdt, indicative of the largest 
inflow rate per unit length, is reached inside the UHR of the oval, 
at r~100 pc. The integrated inflow rate measured at r~50 pc, 
the position closest to the AGN where dM/dt can be estimated 
down to the limit imposed by the spatial resolution of our obser- 
vations, is dM/dt — 8xl0~ 3 M© yr -1 . This inflow rate can sat- 
isfactorily explain the AGN activity of NGC 4579: the reported 
value of dM/dt lies within the range of the mass accretion rates 
derived from the typical bolometric luminosities in low luminos- 
ity AGNs (~10 -2 to ~ 10~ 5 M G yr" 1 from Seyfert to LINERs; see 
compilation by Jogee l200"6b . 



9. Summary and conclusions 

In this paper we have studied the efficiency of angular momen- 
tum transport in the LINER/Seyfert 1 .9 spiral NGC 4579 using a 
complete gravity torque map of the disk of the galaxy. Gas flow 
rates from the outer disk region at r ~ 15 kpc down to the inner 
disk at ~a few 10 pc are derived to search for signatures of secu- 
lar evolution in the fueling process. The CO maps obtained with 
the PdBI and the 30m telescope as part of the NUGA survey are 
used to derive the gravity torque budget in the circumnuclear re- 
gions. Gravity torques on the outer disk gas are derived using the 
HI map obtained with the VLA. The stellar potential is derived 
from a ^-band wide field image of the galaxy. We summarize 
the main results of our study as follows: 

- By analyzing the distribution and kinematics of the gas and 
their relation to other tracers of the ISM and the stellar struc- 
ture we identify signs of dynamical decoupling of several 
patterns on different spatial scales of the disk. The ^T-band 
image reveals a stellar bar of diameter D~12 kpc and mod- 
erate strength, as well as a nested weak oval distortion in the 
inner r -200 pc of the disk. The bar and the oval are not 
aligned, a signature of dynamical decoupling. Molecular gas 
in the inner r <2 kpc disk is distributed in two spiral arms, 
an outer ring, and a central lopsided disk. The morphology 
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of the outer disk revealed by HI shows that the neutral gas is 
currently piling up in a pseudo-ring at r >4 kpc. The pseudo- 
ring is formed by two winding spiral arms that are morpho- 
logically decoupled from the bar structure; this suggests that 
the spirals and the bar do not share a common pattern speed. 

- We have estimated the pattern speeds and principal reso- 
nances of the bar and the oval in the disk of the galaxy. In 
the scenario of dynamical decoupling, the oval has a higher 
pattern speed compared to the primary bar. The gas response 
in the circumnuclear regions traced by the CO maps fol- 
lows the expected gas flow pattern induced by the bar po- 
tential in the presence of two ILRs at r -500 pc (ilLR) and 
r ~ 1 .3 kpc (oILR). We can explain the decoupling of the oval 
with respect to the large-scale bar, as the corotation of the 
oval would be close to the oILR of the large-scale bar. The 
pseudo-ring detected in HI is located inside the bar corota- 
tion (rcR ~6 kpc) and close to the predicted position of the 
UHR of the bar (r UHR -3.8 kpc). 

- The derived gravity torque budget in NGC 4579 shows that 
the fueling process is at work on different spatial scales in 
the disk. In the outer disk, the decoupling of the spiral allows 
the gas to efficiently populate the UHR region inside corota- 
tion and thus produce net gas inflow on intermediate scales. 
This suggests that the corotation barrier of the bar has been 
overcome due to secular evolution processes. The gas stops 
close to the UHR of the bar, where torques become negligi- 
ble. The gas in the inner disk is efficiently funneled by grav- 
ity torques from the oILR down to r -300 pc. Closer to the 
AGN (r<200 pc), gas feels negative torques due to the com- 
bined action of the large-scale bar and the inner oval. The 
two m = 2 modes act in concert to produce net gas inflow 
down to the spatial resolution of our observations, providing 
a clear smoking gun evidence of fueling. 

- In the inner disk, the mass feeding budget is systemati- 
cally negative from r~2 kpc down to r~400 pc. The in- 
flow rate integrated on these scales ( — 1.5 M© yr" 1 ) is 
an order of magnitude higher than the gas feeding rate re- 
quired to sustain the current SFR in this region. Closer to 
the AGN (r<200 pc), the torques are negative down to 
r ~50 pc. The integrated inflow rate measured at r -50 pc 
is — 8x1 0~ 3 M© yr -1 . This inflow rate can satisfactorily ex- 
plain the AGN activity of NGC 4579. 

NGC 4579 is a representative example of a LLAGN in which 
a hierarchy of mechanisms can act in concert to efficiently trans- 
port the gas from the outer disk (on kpc scales) down to vicin- 
ity of the AGN (on ~a few 10 pc scales) through the action of 
gravity torques. The decoupling of dynamical modes on several 
spatial scales (outer spiral, bar, oval) makes it possible to drain 
the angular momentum of the gas and feed the star formation 
in the disk, and possibly the AGN activity itself. In the last step 
down to the AGN, viscosity can also play a significant role in the 
fueling process (e.g., see discussion in GB05). In particular, we 
estimate that the efficiency of viscosity at a typical distance of 
r ~50 pc around the central engine of NGC 457^3 may be com- 
parable to the efficiency of gravity torques. Both mechanisms 
can combine their efforts and produce sufficient AGN feeding 
on extremely short dynamical time-scales of ~l-2 rotation peri- 
ods. 

Acknowledgements. We acknowledge the IRAM staff from the Plateau de Bure 
and from Grenoble for carrying out the observations and help provided during 



5 originally estimated by GB05 at r ~200 pc 



the data reduction. We thank Roberto Neri for fruitful discussions. I.M. acknowl- 
edges financial support from the Spanish Ministerio de Ciencia y Tecnologfa 
(grant AYA2007-62190) and from the Junta de Andalucfa (grant TIC-144). 



References 

Anderson, J. M., Ulvestad, J. S., & Ho, L. C. 2004, ApJ, 603, 42 
Athanassoula, E. 1992, MNRAS, 259, 345 
Baker, A. J. 2000, Ph.D. Thesis 

Baker, A. J., Jogee, S., Sakamoto, K., & Scoville, N. Z. 2003, Active Galactic 

Nuclei: From Central Engine to Host Galaxy, 290, 479 
Barth, A. J., Reichert, G. A., Filippenko, A. V. et al. 1996, AJ, 112, 1829 
Barth, A. J., Ho, L. C, Filippenko, A. V., Rix, H.-W., & Sargent, W. L. W. 2001, 

ApJ, 546, 205 

Bohlin, R. C, Savage, B. D., & Drake, J. F. 1978, ApJ, 224, 132 

Boker, T., Calzetti, D., Sparks, W. et al. 1999, ApJS, 124, 95 

Boone, E, Baker, A. J., Schinnerer, E. et al. 2007, A&A, 471, 113 

Bournaud, E, & Combes, F. 2002, A&A, 392, 83 

Braine, J., & Combes, F. 1992, A&A, 264, 433 

Buta, R., & Combes, F. 1996, Fund. Cosmic Phys., 17, 95 

Canzian, B. 1993, ApJ, 414, 487 

Cardelli, J. A., Clayton, G. C, & Mathis, J. S. 1989, ApJ, 345, 245 

Casasola, V., Combes, E, Garcfa-Burillo, S. et al. 2007, ArXiv e-prints, 712, 

larXiv:0712.0294l 
Chemin, L., Balkowski, C, Cayatte, V. et al. 2006, MNRAS, 366, 812 
Combes, F. 1988, in Galactic and Extragalactic Star Formation, ed. by R. E. 

Pudritz & M. Fiched, NATO Advanced Science Institutes (ASI) Series C, 

Volume 232, p. 475 

Combes, F. 2003, in Active Galactic Nuclei: from Central Engine to Host Galaxy, 
ed. by S. Collin, F. Combes & I. Shlosman. ASP (Astronomical Society of 
the Pacific), Conference Series, V ol. 290, p. 411. 
Combes, F. 2007, ArXiv e-prints, 709, arXiv:0709.0091 
Combes, E, & Sanders, R. H. 1981, A&A, 96, 164 
Combes, E, Garcfa-Burillo, S., Boone, F. et al. 2004, A&A, 414, 857 
Comeron, S., Knapen, J. H., Beckman, J. E., & Shlosman, I. 2008, A&A, 478, 
403 

Contini, M. 2004, MNRAS, 354, 675 

Daigle, O., Carignan, C, Amram, P. et al. 2006, MNRAS, 367, 469 
Dewangan, G. C, Griffiths, R. E., Di Matteo, T., & Schurch, N. J. 2004, ApJ, 
607, 788 

Eracleous, M., Shields, J. C, Chartas, G., & Moran, E. C. 2002, ApJ, 565, 108 
Eskridge, P. B., Frogel, J. A., Pogge, R. W. et al. 2002, ApJS, 143, 73 
Ferrarese, L., & Ford, H. 2005, Space Sci. Rev., 116, 523 
Filippenko, A. V, & Sargent, W. L. W. 1985, ApJS, 57, 503 
Garcfa-Burillo, S., Guelin, M., & Cernicharo, J. 1993, A&A, 274, 123 
Garcfa-Burillo, S., Sempere, M. J., Combes, E, & Neri, R. 1998, A&A, 333, 864 
Garcfa-Burillo, S., Combes, E, Eckart, A. et al. 2003a, in ASP Conf. Ser., Active 

Galactic Nuclei: from Central Engine to Host Galaxy, ed. by S. Collin, F. 

Combes,& I. Shlosman, 423. 
Garcfa-Burillo, S., Combes, E, Hunt , L. K. et al. 2003b, A&A, 407, 485 
Garcfa-Burillo, S., Combes, E, Schinnerer, E., Boone, E, & Hunt, L. K. 2005, 

A&A, 441, 1011 

Gil de Paz, A.; Boissier, S.; Madore, B. E; Seibert, M.; Joe, Y. H. et al. 2007, 
ApJS, 173, 185 

Guilloteau, S., & Lucas, R. 2000, in ASP Conf. Ser.: Imaging at Radio through 
Submillimeter Wavelengths, ed. by J. G. Mangum & S. J. E. Radford, vol. 
299 

Guhathakurta, P., van Gorkom, J. H., Kotanyi, C. G., & Balkowski, C. 1988, AJ, 
96, 851 

Haan, S., Schinnerer, E., Mundell, C. G., Garcfa-Burillo, S., & Combes, F. 

2008a, AJ, 135, 232 
Haan, S. et al. 2008b, ApJ, submitted 
Ho, L. C, & Ulvestad, J. S. 2001b, ApJS, 133, 77 
Ho, L. C, Feigelson, E. D., Townsley, L. K et al. 2001a, ApJ, 549, L51 
Ho, L. C, Filippenko, A. V, & Sargent, W. L. W. 1997, ApJS, 112, 315 
Holtzman, J. A., Burrows, C. J., Casertano, S., et al. 1995, PASP, 107, 1065 
Hopkins, P. E, & Hernquist, L. 2006, ApJS, 166, 1 

Hummel, E., van der Hulst, J. M., Keel, W. C, & Kennicutt, R. C, Jr. 1987, 
A&AS,70, 517 

Hunt, L. K., Combes, E, Garcfa-Burillo, S. et al. 2008, A&A, 482, 133 
Hutchings, J. B., & NefF, S. G. 1992, AJ, 104, 1 

Jarrett, T. H., Chester, T., Cutri, R., Schneider, S. E., & Huchra, J. P. 2003, AJ, 
125, 525 

Jogee, S. 2006, Physics of Active Galactic Nuclei at all Scales, 693, 143 
Jogee, S., Baker, A. J., Sakamoto, K. et al. 2001, in ASP Conf. Ser.:The Central 
Kiloparsec of Starbursts and AGN: The La Palma Connection, vol 249, 612 
Kennicutt, R. C, Jr. 1998, ARA&A, 36, 189 



S. Garcfa-Burillo et al.: Molecular gas in NUclei of GAlaxies (NUGA) 



21 



Kenney, J. D. P., & Young, J. S. 1989, ApJ, 344, 171 
Kennicutt, R. C, Jr., Armus, L., Bendo, G. et al. 2003, PASP, 115, 928 
King, A. R., & Pringle, J. E. 2007, MNRAS, 377, L25 
Knapen, J. H., Shlosman, L, & Peletier, R. F. 2000, ApJ, 529, 93 
Knapen, J. H., de Jong, R. S., Stedman, S., & Bramich, D. M. 2003, MNRAS, 
344, 527 

Kohno, K., Matsushita, S., Vila-Vilaro, B. et al. 2001 in ASP Conf. Ser.:The 
Central Kiloparsec of Starbursts and AGN: The La Palma Connection, vol 
249, 672 

Koopmann, R. A., Kenney, J. D. P., & Young, J. 2001, ApJS, 135, 125 

Krause, M., Fendt, C, & Neininger, N. 2007, A&A, 467, 1037 

Krips, M., Eckart, A., Neri, R. et al. 2005, A&A, 442, 479 

Krips, M., Eckart, A., Krichbaum, T. P. et al. 2007, A&A, 464, 553 

Kuno, N., Sato, N., Nakanishi, H. et al. 2007, PASJ, 59, 117 

Lindt-Krieg, E., Eckart, A., Neri, R. et al. 2008, A&A, 479, 377 

Maoz, D., Koratkar, A., Shields, J. C. et al. 1998, AJ, 116, 55 

Maoz, D., Nagar, N. M., Falcke, H., Wilson, A. S. 2005, ApJ, 625, 699 

Marquez, I., Durret, E, Masegosa, J. et al. 2000, A&A, 360, 431 

Martini, P. 2004 in Coevolution of Black Holes and Galaxies, ed. by L. C. Ho, 

Cambridge University Press, p. 170. 
Mulchaey, J. S., & Regan, M. W. 1997, ApJ, 482, L135 

Pahre, M. A., Ashby, M. L. N., Fazio, G. G, & Willner, S. P. 2004, ApJS, 154, 
235 

Peeples, M. S., & Martini, P. 2006, ApJ, 652, 1097 
Pogge, R. W. 1989, ApJS, 71, 433 

Pogge, R. W., Maoz, D., Ho, L. C, & Eracleous, M. 2000, ApJ, 532, 323 

Regan, M. W., & Teuben, P. 2003, ApJ, 582, 723 

Regan, M. W., & Teuben, P. J. 2004, ApJ, 600, 595 

Rubin, V. C, Waterman, A. H., & Kenney, J. D. P. 1999, AJ, 118, 236 

Sellwood, J. A., & Sparke, L. S. 1988, MNRAS, 231, 25P 

Sanders, D. B., Mazzarella, J. M., Kim, D. -C, Surace, J. A., & Soifer, B. T. 

2003, AJ, 126, 1607 
Solomon, P. M., & Barrett, J. W. 1991 in Dynamics of Galaxies and their 

Molecular Cloud Distributions, from the IAU Symp. 146, ed by F. Combes 

& F. Casoli, Kluwer Academic Publishers, 235. 
Stauffer, J. R. 1982, ApJ, 262, 66 

Terashima, Y, Ho, L. C, Ptak, A. F. et al. 2000, ApJ, 535, L79 

Tully, R. B., & Fisher, J. R. 1988, Catalog of Nearby Galaxies, by R. Brent 

Tully and J. Richard Fisher, pp. 224. Cambridge, UK: Cambridge University 

Press, April 1988 
Ulvestad, J. S., & Ho, L. C. 2001, ApJ, 562, L133 

Wada, K. 2004 in Coevolution of Black Holes and Galaxies, ed. by L. C. Ho, 
Cambridge University Press, p. 187 



